Modeling Absolute Redox Potentials of Ferrocene in the Condensed Phase.

Absolute thermodynamic quantities for critical chemical reactions are needed to determine the role of solvents and reactive environments in catalysis and electrocatalysis. Theoretical methods can provide such quantification but are often hindered by the innate complexity of electron correlation and dynamic relaxation of solvent environments. We present and validate a protocol for calculating the redox potentials of the ferrocene/ferrocenium redox pair in acetonitrile. Equation-of-motion and effective fragment potential (EFP) methods are used to characterize the adiabatic and vertical ionization potentials as well as the electron affinity processes. We benchmark molecular mechanics against the EFP model to show the differences in the ferrocene electronic polarizability in two redox states. Our best estimate of the redox potential (4.94 eV) agrees well with the experimental value (4.93 eV). This demonstrates the ability of modern computational methods to predict absolute redox potentials quantitatively and to quantify the correlation of dynamic effects, which underlie their origin.

Binding and stability of MgO monomers on anatase TiO2(101).

In catalysis, MgO is often used to modify the acid-base properties of support oxides and to stabilize supported metal atoms and particles on oxides. In this study, we show how the sublimation of MgO powder can be used to deposit MgO monomers, hither on anatase TiO2(101). A combination of x-ray electron spectroscopy, high-resolution scanning tunneling microscopy, and density functional theory is employed to gain insight into the MgO monomer binding, electronic and vibrational properties, and thermal stability. In the most stable configuration, the Mg and O of the MgO monomer bind to two surface oxygens and one undercoordinated surface titanium, respectively. The additional binding weakens the Mg-O monomer bond and makes Mg more ionic. The monomers are thermally stable up to 600 K, where the onset of diffusion into the TiO2 bulk is observed. The monomeric MgO species on TiO2(101) represent an ideal atomically precise system with modified acid-base properties and will be employed in our future catalytic studies.

Creating self-assembled arrays of mono-oxo (MoO3)1 species on TiO2(101) via deposition and decomposition of (MoO3)n oligomers.

Hierarchically ordered oxides are of critical importance in material science and catalysis. Unfortunately, the design and synthesis of such systems remains a key challenge to realizing their potential. In this study, we demonstrate how the deposition of small oligomeric (MoO3)1-6 clusters-formed by the facile sublimation of MoO3 powders-leads to the self-assembly of locally ordered arrays of immobilized mono-oxo (MoO3)1 species on anatase TiO2(101). Using both high-resolution imaging and theoretical calculations, we reveal the dynamic behavior of the oligomers as they spontaneously decompose at room temperature, with the TiO2 surface acting as a template for the growth of this hierarchically structured oxide. Transient mobility of the oligomers on both bare and (MoO3)1-covered TiO2(101) areas is identified as key to the formation of a complete (MoO3)1 overlayer with a saturation coverage of one (MoO3)1 per two undercoordinated surface Ti sites. Simulations reveal a dynamic coupling of the reaction steps to the TiO2 lattice fluctuations, the absence of which kinetically prevents decomposition. Further experimental and theoretical characterizations demonstrate that (MoO3)1 within this material are thermally stable up to 500 K and remain chemically identical with a single empty gap state produced within the TiO2 band structure. Finally, we see that the constituent (MoO3)1 of this material show no proclivity for step and defect sites, suggesting they can reliably be grown on the (101) facet of TiO2 nanoparticles without compromising their chemistry.

Polarizable embedding for simulating redox potentials of biomolecules.

Redox reactions play a key role in various biological processes, including photosynthesis and respiration. Quantitative and predictive computational characterization of redox events is therefore highly desirable for enriching our knowledge on mechanistic features of biological redox-active macromolecules. Here, we present a computational protocol exploiting polarizable embedding hybrid quantum-classical approach and resulting in accurate estimates of redox potentials of biological macromolecules. A special attention is paid to fundamental aspects of the theoretical description such as the effects of environment polarization and of the long-range electrostatic interactions on the computed energetic parameters. Environment (protein and the solvent) polarization is shown to be crucial for accurate estimates of the redox potential: hybrid quantum-classical results with and without account for environment polarization differ by 1.4 V. Long-range electrostatic interactions are shown to contribute significantly to the computed redox potential value even at the distances far beyond the protein outer surface. The approach is tested on simulating reduction potential of cryptochrome 1 protein from Arabidopsis thaliana. The theoretical estimate (0.07 V) of the midpoint reduction potential is in good agreement with available experimental data (-0.15 V).

Extension of the Effective Fragment Potential Method to Macromolecules.

The effective fragment potential (EFP) approach, which can be described as a nonempirical polarizable force field, affords an accurate first-principles treatment of noncovalent interactions in extended systems. EFP can also describe the effect of the environment on the electronic properties (e.g., electronic excitation energies and ionization and electron-attachment energies) of a subsystem via the QM/EFP (quantum mechanics/EFP) polarizable embedding scheme. The original formulation of the method assumes that the system can be separated, without breaking covalent bonds, into closed-shell fragments, such as solvent and solute molecules. Here, we present an extension of the EFP method to macromolecules (mEFP). Several schemes for breaking a large molecule into small fragments described by EFP are presented and benchmarked. We focus on the electronic properties of molecules embedded into a protein environment and consider ionization, electron-attachment, and excitation energies (single-point calculations only). The model systems include chromophores of green and red fluorescent proteins surrounded by several nearby amino acid residues and phenolate bound to the T4 lysozyme. All mEFP schemes show robust performance and accurately reproduce the reference full QM calculations. For further applications of mEFP, we recommend either the scheme in which the peptide is cut along the Cα-C bond, giving rise to one fragment per amino acid, or the scheme with two cuts per amino acid, along the Cα-C and Cα-N bonds. While using these fragmentation schemes, the errors in solvatochromic shifts in electronic energy differences (excitation, ionization, electron detachment, or electron-attachment) do not exceed 0.1 eV. The largest error of QM/mEFP against QM/EFP (no fragmentation of the EFP part) is 0.06 eV (in most cases, the errors are 0.01-0.02 eV). The errors in the QM/molecular mechanics calculations with standard point charges can be as large as 0.3 eV.